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We investigate the charge transport in one-dimensional arrays of Josephson junctions. In the 
interesting regime of "small charge solitons" (polarons), AEj > Ec > Ej, where A is the (elec- 
trostatic) screening length, the charge dynamics is strongly influenced by the polaronic effects, i.e., 
by dressing of a Cooper pair by charge dipoles. In particular, the soliton's mass in this regime 
scales approximately as EJ 2 . We employ two theoretical techniques: the many body tight-binding 
approach and the mean- field approach. Results of the two approaches agree in the regime of "small 
charge solitons" . 



I. INTRODUCTION 



Physics of one- and two-dimensional arrays of Joseph- 
son junctions is surprisingly rich. Both 1-D and 
2-D 10-14 arrays (including granulated superconducting 
films) have been extensively investigated. Yet, many 
unanswered questions remain. In particular, transport 
properties of 1-D arrays of Josephson junctions are still 
not fully understood. Experiments 4,7,8 show various phe- 
nomena related to superconductor-insulator transitions, 
Coulomb blockade, hysteresis, mixed Josephson-quasi- 
particle effects etc.. One of the challenging questions is 
the value and the origin of the mass of the charge carriers 
in the insulating regime. In the theoretical studies of Her- 
mon et al. it was shown that, if the grains have a large 
kinetic (or geometric) inductance, the system's dynam- 
ics are governed by the sine-Gordon model and, therefore, 
kink-like topological excitations, i.e., charge solitons, are 
the charge carriers. In Ref. the domain of applicability 
of this sine-Gordon description was analyzed. Simultane- 
ous experiments by Haviland and Delsing demonstrated 
the Coulomb blockade in 1-D arrays of JJs consistent 
with the existence of charge solitons. In the later exper- 
iments of Haviland's group ' considerable hysteresis in 
the I-V characteristic of the array was observed and at- 
tributed to a very large kinetic inductance. The physical 
origin of this inductance remained unclear. A few years 
later, Zorin 15 pointed out that a current biased small- 
capacitance JJ develops an inductive response on top of 
the capacitive one. This phenomenon was called Block 
inductance. A closely related inductive coupling between 
two charge qubits was studied in Ref. . The role of 
the Bloch inductance in Josephson arrays was studied in 



Ref. for the case of an infinite screening length, i.e., 
when the array serves as a zero-dimensional lumped cir- 
cuit element. 

In this paper we employ two complimentary techniques 
to study the charge propagation in infinite Josephson ar- 
rays with finite but large screening length. We consider 
arrays free of disorder. Specifically we concentrate on cal- 
culating the effective mass of the charge carriers. Both 
approaches, the many-body tight-binding technique and 
the mean-field technique agree for not very small ratios 
Ej/Ec- In particular, the effective mass of a charge soli- 
ton scales approximately as EJ 2 in this regime. For full 
transport description one has to treat the effects of the 
array's boundaries as well as those of the disorder. Yet, 
our result about the effective mass is clearly relevant for 
further investigation of the transport. 

II. THE SYSTEM 

The system under study is shown in Fig. 1. The 
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FIG. 1: Array of Josephson junctions. 



Josephson junctions with capacitance C connect the su- 
perconducting grains to each other and each grain has a 
capacitance Co to the ground. Typical values are C ~ 
IfF and Cq ~ 5 — 20 aF. The system is governed by 
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the usual Hamiltonian consisting of the Coulomb charg- 
ing energy (kinetic energy) and the Josephson tunneling 
(potential energy): 

H= l -Y J U(r-r , )n r n r/ -EjY J ^(O r -e r _ 1 ) . (1) 



Here n r are integer-valued island charges (in units of 
2e) and r are the corresponding canonically conjugate 
phases, [n r ,e^'] = e %er 8 rr >. The matrix of Coulomb 
interaction U (r) is given by 



U(r) = 2E C 



f 

J — 7T 



dk 



a ikr 



2ttA- 2 -2(cosA;-1) ' 



(2) 



with Ec = (2e 2 )/2C being the charging energy and 
A = ^/CjCo the screening length which determines the 
spatial extent of the Coulomb interaction. In this paper 
we consider A ^> 1. 



III. TIGHT-BINDING APPROACH 
A. Qualitative discussion 

In this section we explore the properties of Josephson 
arrays in the Coulomb blockade regime Ej <C Ec- We 
consider the sector of the Hilbert space with exactly one 
extra Cooper pair in the array. The simplest (and hav- 
ing minimal charging energy) representative of the unit 
charge sector is the state in which the extra Cooper pair 
resides on some island R with all the other islands be- 
ing neutral. The charging energy of such a state is given 
by /io = \U{G) ~ AEc/2. This is approximately the 
energy (rather high!) one has to invest to insert one 
Cooper pair into the array. Once the Cooper pair has 
been inserted it is free to move from one site to its neigh- 
bor via the hopping provided by the Josephson part of 
the Hamiltonian. In the limit of vanishingly small Ej 
only the simplest charge configurations described above 
are important and we are led to the trivial tight-binding 
band E(k) = — Ej cos k for an extra Cooper pair in the 
Josephson chain (cf. 2 ' 3 ). 

The peculiarity of the ID- Josephson chain, first no- 
ticed in Refs. 2 ' 3 and used in Ref. , is that the sim- 
ple picture sketched above is valid only for extremely 
small Ej < Ec/ 'A. The reason is the presence of a 
large number of states lying at small energy ~ Ec / 'A 
above the basic states (as opposed to much larger en- 
ergy Ec which one might expect and which indeed hap- 
pens in higher dimensions). One particular example 
is the charge configuration |1,— 1,1) (Cooper pair and 
a properly oriented dipole nearby) having the energy 
|Z7(0) - 2/7(1) + U{2) « /i + E c /A. Thus, in the pa- 
rameter range AEj > Ec > Ej called by the authors 
of Ref. the small soliton regime, Cooper pair inserted 
into the chain gets strongly dressed by virtual dipoles and 
the simplest tight-binding scheme breaks down. Dipole 



dressing was also mentioned in the context of transport 
in ion channels 18 . 

In reference the properties of the small charge soli- 
tons were addressed by successive inclusion of the charge 
configurations (up to 32 states) with larger and larger 
energies into the tight-binding scheme. A similar scheme 
was developed for polarons in Ref. . In this paper we 
construct a comprehensive description of the low-lying 
(with energies much smaller Ec) states in terms of a 
particular spin-1/2 model. We derive an effective Hamil- 
tonian governing the model dynamics within the low en- 
ergy subspace. We then develop a tight-binding approach 
with arbitrary number of the charge states taken into ac- 
count. 



B. Structure of the low energy subspace 

Let us first define more precisely what we mean under 
the low lying states in the sector with total charge 1 
and construct the complete classification of these states. 
Let us consider some charge configuration of size w. It is 
clear that the energy of such a configuration will certainly 
exceed Ec if w > A (from now on we count energies 
from the energy /i of a single Cooper pair). Thus for 
the low-lying configurations w < A and we can expand 
the charging energy in powers of 1/ A as 



H c = -^tt y^J r ~ r r \n r n r t + ^ y^(r - r f ) 2 n r n r > + 



E c 



4A 



(3) 

We see that the typical charge configurations have large 
energy ~ Ec ^> Ej and are not important for the low 
energy physics. The exceptions are the states nullifying 
the first term in Eq. (3) and having the energy ~ wEc / A. 
Note, that the first term of the expansion (3) can not take 
negative values. Otherwise there would exist configura- 
tions with electrostatic energy smaller than /io- As long 
as w < AEj/Ec these charge configurations hybridize 
effectively with the basic one leading to the formation of 
the small charge soliton. Thus, the condition 



\r — r'\n r n r > 







(4) 



is the mathematical definition of the low energy sub- 
space in the unit charge sector (we call it also the proper 
space) . It can be shown (Appendix A) that the subspace 
(4) consists of all the configurations with two properties: 
a) all the islands' charges n r equal ±1 or 0; b) any two 
charged islands separated by an arbitrary number of neu- 
tral islands have opposite charges. For example, the con- 
figurations 1 1,0,— 1,1) and 1 1, 0, — 1,0, 1) belong to the 
low-energy space while the configurations |2,0, — 1) and 
1 1, 0, 1,0, — 1) do not. To describe the proper subspace in 
a more clear way let us introduce variables a r defined on 
the links of the chain (the link r is the link connecting 
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FIG. 2: Connection between <j z r (denned on the links) and the 
charges of the islands n r . In the low-energy subspace a* can 
take only two values ±1 and are ^-projections of a spin- 1/2. a) 
The basic charge configuration with only one charged island 
corresponds to an abrupt domain wall in terms of a* . b) More 
complicated charge configuration (two additional dipoles) cor- 
responding to a domain wall of finite thickness in the spin 
language. 



islands r and r + 1) 



./> r _l_i 



r' <r 



(5) 



The connection between variables a* and charges n r 
for two configurations in the low-energy subspace is illus- 
trated in Fig. 2. From the definition of o z r and the prop- 
erties of the states in the low-energy subspace one imme- 
diately concludes that the low-energy configurations are 
described by o z r = ±1 for all r, i.e., the low-energy sub- 
space is isomorphic to the space of states for a spin- 1/2 
chain with a* being the z-projections of the spins. 

Due to the constraint n r = 1 the variable a* satis- 
fies the boundary conditions 



-1 



>+oo 



1 . 



(6) 



Thus, extra Cooper pair in the chain is described by a 
domain wall in the spin language. 



C. Projecting the Hamiltonian 

Having understood the structure of the low-energy 
space of the model we can project the full Hamiltonian 
(1) onto the proper subspace. The projection is carried 
out by noting that Cooper pair tunneling between two 
neighboring islands corresponds to the spin flip in the link 
between them. Thus the Josephson part of the Hamilto- 
nian is given by 



Hj -Ej ^2 a r 



(7) 



Rewriting the charging energy in terms of spin variables 
we arrive at 



H 



1 



-i)U(r- 



(8) 

To determine the spectrum of the single-charge sector of 
the Hamiltonian (8) we impose the boundary conditions 
(6) indicating the presence of a domain wall. 

The Hamiltonian (8) takes into account the low en- 
ergy charge configurations of arbitrary width w. We 
understand however that the configurations with w ^> 
AEj/Ec are not important at low energies. Thus, we can 
further reduce the phase space by dropping out all the 
configurations of the width w larger than some wq. We 
expect that at wo ^> AEj/Ec the resulting low energy 
states are independent of wo and approximate correctly 
those of Hamiltonian (8). 

Any state containing a domain wall of the width less 
than wq is completely specified by the position R of 
the first spin up (which we call the coordinate of the 
charge soliton or domain wall) and the values of the z- 
projections of the next wq spins {ai, . . . cr Wo }. Given the 
state 



|i?)|ai,...<j Wo ) 



(9) 



one can reconstruct the z-projections of all spins in the 
chain according to 



-1, 
_ 1, 

1, 



R 



r < R , 
r = R , 
- 1 < r < R + w , 
r > R + iu . 



(10) 



For example, if we choose wo — 5 the states shown on 
Fig. 2a) and 2b) can be written as 



|iJ = 0)|t,t,t,t,t) , 

\R=-2)\i,i,tAA) ■ 



(11) 

(12) 



In Appendix B we describe how to project the Hamil- 
tonian (8) onto the space of configurations with sizes less 
or equal than i^o- We also perform a transition from the 
coordinates \R) to the quasi- momentum k. The result 
reads 



H = H c -EjY J °r 

r=l 

w w 



-E 



r=l 



r=l 



(13) 



where T is the operator of the right cyclic shift defined 
by T\ai , . . .a WQ ) = \a Wo , a t , . . ,cr 2 ). 

The Hamiltonian (13) constitutes the main result of 
this section. For wq — 1, 2, ... 5 it can be shown to pro- 
duce results equivalent to that of reference 17 . Equation 
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(13) reduces the initial many-body problem to a finite di- 
mensional Hamiltonian, readily accessible to numerics as 
long as not too large (wq < 20) charge configurations are 
important. In the next sections we present the results of 
numerical analysis of the Hamiltonian (13) and compare 
the results to those of the mean-field approach. 



D. Results of the tight-binding approach 

An example of the band structure obtained within the 
tight-binding approach is shown in Fig. 3. In Fig. 4 
the two lowest bands are shown. We observe that the 



E'(k)/£, 




FIG. 3: Spectrum of the tight-binding Hamiltonian (13) for 
Ej/Ec = 0.4 and A = 10. The number of charge states taken 
into account equals 2 7 . 




FIG. 5: Solid (blue) line: the group velocity corresponding 
to the lowest energy band of Fig 4. Dashed (red) line: the 
naive tight-binding group velocity with no dressing by dipoles 
taken into account. 



of the charge carriers. In the tight-binding approach we 
define raxB = ^ 2 ( ~ ) , where E (k) is the 



dispersion of the lowest band (ground state). In what 
follows we will compare this mass with the results of the 
mean-field theory. 



1. Persistent current 

As a first obvious application of our results consider a 
ring-shaped array of N junctions with exactly one extra 
Cooper pair in it. If an external magnetic flux <& ex t is 
applied a persistent current will emerge. The periodic 
boundary condition for the Bloch wave with wave vector 
k reads 



mvE; 




FIG. 4: Two lowest energy bands for A = 10 and Ej/Ec — 
0.4. The number of charge states taken into account equals 



lowest band is parabolic for small momenta k and flat- 
tens in the outer part of the Brillouin zone. This phe- 
nomenon was already observed in Ref. 17. To further 
emphasize the dispersion relation of the lowest band in 
Fig. 5 we show the group velocity of the soliton (dressed 
Cooper pair) as compared to the one of an undressed 
Cooper pair. We find that the flattening of the disper- 
sion relation in the outer region of the Brillouin zone 
leads to zero group velocity. In this paper we concen- 
trate mostly on the investigation of the effective mass 



AkN 



2tt^ 



(14) 



Thus, as the external flux varies between — <&q/2 and 
<E>o/2, the relevant wave vector varies between —tt/N 
and tt/N. For large enough N the interval [—tt/N, tt/N] 
is safely within the domain of parabolic dispersion re- 
lation. Thus we use the effective mass approximation 
and obtain for the persistent current in the interval 

$ext e [-$ /2,$o/2] 



2e hk 



N m, 



2e 2nh§ e 



eff 



N 2 m 



eff 



$0 



(15) 



where m e f f is the effective mass of the charge carrier (in 
the tight-binding approach we obtained m e ff = tutb)- 
Thus, the amplitude of the persistent current oscillations 
is given by 



2irhe 



N 2 m, 



eff 



(16) 



With no polaronic effects taken into account, i.e., for a 
bare Cooper pair we would have E b Q are {k) = -Ej cos k 
h 2 /Ej. Thus we obtain 



and m^Yf 



Io 



2neEj m b e ff 



Y ^bare 
TTic m eff 



HN 2 



m 



eff 



N 2 m, 



eff 



(17) 
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where I c is the critical current of a single Josephson junc- 
tion. We observe that the effective mass reduction via the 
polaronic effects enhances the persistent current. 



IV. MEAN-FIELD THEORY 

A. Description in terms of continuous polarization 
charges 

An alternative description of the charge propagation 
in the array is given in terms of the continuous po- 
larization charges, e.g., the screening charges q& ate on 
the gate capacitances Co (see Fig. 6). For the system 
described in the previous section the continuous polar- 
ization charges are enslaved to the discrete charges n r . 
That is, once a tunneling process occurs and the distri- 
bution n r changes, the polarization charges adjust im- 
mediately to the new situation. To allow formally in- 
dependent dynamics of polarization charges we intro- 
duce infinitesimal inductances Lq as shown in Fig. (6). 
This leads to two independent degrees of freedom per 



L 
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FIG. 6: Array of Josephson junctions with infinitesimal in- 
ductances Lq. 

cell of the array. One quantized charge degree of free- 
dom m r = XlfcLr n k is the number of Cooper pairs that 
have tunneled through junction number r. Its conjugate 
phase is given by <\> T = 6 r — 6 r -i and the commutation 
relations read [m^e^'j = e % ^ r 8 rr >. The second contin- 
uous charge degree of freedom Q r = ^ q g r ? te + 2em_oo 

r' <r 

is equal to the polarization charge that has arrived at 
the junction number r or, alternatively, the integral of 
the displacement current flowing into junction r. The 
conjugate variable <l> r is the magnetic flux on inductance 
Lq in cell number r. The commutation relation reads 
[§ ri Q r >] = ih5 rr > . We obtain the following Hamiltonian 
of the array 



H = 



(2em r — Q r ) 2 



2C 

(Qr — Qr-l) 2 

2C 



Ej cos(^ r ) 



$2 

^ r 

2L 



(18) 



B. Mean-field approximation 

The mean-field description is based on the Heisenberg 
equations of motion for the polarization charge Q r fol- 



lowing from (18): 



Ln 



C 



(Q r — 2em r ) 



2Qr Qr — 1 Qr-\-\ 

Co 



(19) 



We average Eq. (19) over the state of the system and 
obtain 



L (Q r ) = -V r 



2{Q r ) ~ (gr+l) ~ (Qr-l) 

Co 



(20) 



where V r = (-^(Qr — 2em r )) is the expectation value 
of the voltage drop across junction number r. In the 
mean-field approximation we calculate V r by replacing 
the operators Q r by their average values (Q r )(t) in the 
Hamiltonian (18). Thus the problem factorizes to many 
single-junction ones. Each junction is governed by the 
Hamiltonian 



H(Q(t)) 



(2em - Q(t)) 2 
2C 



Ej cos ( 



(21) 



where we have dropped the index r. The gate charge 
Q(t) is a given function of time (to be replaced in each 
junction by (Q r )(t)). For the expectation value of the 
voltage we then obtain V = (OqH). The problem is now 
to find the quantum state of the junction in which the 
average (dqH) should be evaluated. We do so assuming 
that (Q r )(t) is a slow function of time. This assumption 
should be checked for self-consistency later. 

The Hamiltonian (21) possesses the (adiabatic) spec- 
trum with discrete eigenvectors \e n (Q(t))) obeying 
(e n |e m ) = Snm and eigenvalues E n (Q(t)), cf. Fig. 7. 
The general wave function is a superposition \Sfr(t)) = 
^2a n (t) \e n (Q(t))}. Our aim is to determine \^(t)} for 




FIG. 7: The ground and first excited state E , E 1 of (21) 

a given function Q(t). We restrict ourselves to the adia- 
batic case, i.e., we keep only terms of order Q, Q 2 . After 
a calculation presented in Appendix C we arrive at 

V = (d Q H) = d Q E + L B Q + \(d Q L B )Q 2 . (22) 

Here we defined the Bloch inductance first introduced by 

Zorin 15 : 



j o fc 2T (en\d Q H\e ) 

Lb = 2H ^ (E n -E r 



n>0 



(23) 
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(In Ref. only the first excited state (n = 1) in (23) 
was taken into account and the contribution ex 3qLb in 
(22) was omitted.) For Ej <C Eq the Bloch inductance 
Lb is sharply peaked around Q = e (see Fig. 8). In 
the opposite case, Ej ^> Ec, the Bloch inductance Lb is 
nearly constant, Lb ~ Lj = 4 2 ^ . Combining Eq. (22) 



Bloch inductance having a peak value L max = Ls(e) 
at Q — e we obtain the minimal light velocity c m i n = 
/T 1 . ii) The ground state energy is still a 2e 
periodic function of Q but it is no longer proportional 
to cosQ. hi) Since Lb depends strongly on Q, the third 
term of (25) is very important. 




FIG. 8: Solid line: the Bloch inductance L B (Q) for Ej/E c = 
0.5, measured in units of the Josephson inductance Lj. 
Dashed line: Lb(Q) = Lj for Ej/Ec — oo. 

and the self-consistency equation (20) we obtain 



L 



OVr 



^Qr Qr — 1 Qr-\-l 



Co 



d Q E (Q r ) + L B (Q r )Qr + ^d Q L B (Q r )Q 2 r 



(24) 



where we substituted (Q) —> Q for clarity. We observe 
that the kinetic inductance is superseded by the Bloch 
inductance at least around Q = e and we can safely as- 
sume Lq — >• 0. Yet, in the regime Ej <C Eq-, when Lb(Q) 
is exponentially small in the regions Q ~ and Q ~ 2e, a 
finite geometric or kinetic inductance L could be impor- 
tant. In the continuum limit, i.e., after the substitution 
r _i — 2Q r — >> d^Q = Q" ', equation (24) reads 



ir+l 



L B (Q)Q 



1 



Co 



(SqL b )Q 2 



dEp(Q) 
dQ 



= 



(25) 



We now make the very important observation that (25) 
is the equation of motion for the following Lagrangian 
density 



1 



C{Q,Q) = -L B (Q)Q 



1 

2C~o 



Q' 2 -E (Q) 



(26) 



In the limit Ej ^> Ec, when Lb ~ const and Eo(Q) oc 
cos Q we obtain the usual sine-Gordon equation. On the 
other hand, in the limit Ej <C Ec equation (25) differs 
in several aspects from the sine-Gordon equation: i) The 
first two terms of (25) describe a wave guide with a in- 
dependent "light velocity" c(Q) = , 1 With the 



y/L B (Q)C ' 



C. Solitonic solutions 

We are now searching for a solitary wave traveling with 
velocity v by plugging the ansatz Q(r — vt) into Eq. (25). 
This gives the following differential equation 



_9 
dr 



1 



L B (Q)v 2 



1 

Co 



Q ,Z + E (Q) 



= . 



(27) 



Integrating we obtain 



r = ± 



Here, E m ^ n 



Q(r) 

I 

Q(ro) 



dQ 



2Co(E (Q) — E min ) 
l-L B (Q)C v 2 



-1/2 



(28) 



is an integration constant and ± stands for 
the soliton / antisoliton solution. We impose the bound- 
ary conditions Q(— oo) = and Q(+oo) = 2e to describe 
the propagation of a single Cooper pair in the array. 
This also fixes the integration constant, E min = Eq(0) = 
Eo(2e). The solitonic solutions only exist for v < c m i n . 



D. Lorentz contraction 

In the limit Ej ^> Ec, Eq. (25) reduces to the sine- 
Gordon equation and is Lorentz invariant. Thus solitons 
undergo the usual Lorentz contraction. In the other limit, 
Ej <C Ec, Eq. (25) is not Lorentz invariant. The Lorentz 
contraction of the soliton takes a very peculiar shape. 
Consider a soliton moving with velocity v approaching 
Cmin (we postpone for a moment a discussion on whether 
this is consistent with adiabaticity ) . For the center of the 
soliton, where Q ~ e, the relativist ic regime is reached 
and it is Lorentz contracted (see Fig. 9). In contrast, the 
soliton's tales, where Q ~ or Q ~ 2e are unaffected by 
Lorentz contraction. 



E. Rest energy and dynamical mass of the soliton 

Using the Lagrangian density (26) we find the energy 
of a soliton 

2e 



o v 

For small velocities we expand and obtain E so i(v) = 
Erest + \ mkinV 2 + 0(v 4 ). The rest energy of the soliton 
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is given by 



2e 



E r 



dQ\ 



I2[E (Q)-E (0)} 



o 



Co 



For the kinetic mass we obtain 

2e 




(30) 



dQ L B (Q) V2[E (Q) - E o (0)]C o . (31) 



FIG. 10: Solid (red) curve: mean- field mass rrikin- Lower 
dashed (blue) curve: tight-binding mass ttitb- Upper dashed 
(blue) curve: upper boundary for ttitb- 

for Ej/Eq > 0.3. This is the main result of this paper. 
The mass scales approximately as ~ Ej 2 . As the conver- 
gence of the tight-binding approach gets worse with rais- 
ing ratio we also show the uncertainty of the result 
by giving the upper boundary for ttitb (upper dashed 
curve in Fig. 10). 



In the limit Ej ^> Ec-, when Lb ~ Lj = 
const., we obtain as expected the relativist ic relation 
E rest w m fcin c^ in (in this limit the light velocity is Q- 
independent, c(Q) « c min w (LjC )~ 1/2 )- 

In the opposite charging limit, Ej <C Ec, no such 
relativistic relation exists. A simple estimate then gives 



E r 




Ej 



Er 



(32) 



consistent with the result of Sec. Ill A. 

In the limit Ej <C Ec the Bloch inductance Lb(Q) 
is sharply peaked around Q = e and the integral of 
Eq. (31) is dominated by a small vicinity of this point. 
Here a two-state approximation is valid which gives 

L B = 2 (fa) L j sin 5 (0) with cot6> = (^) fa. As 

Lb{Q) is sharply peaked around Q = e, we can replace 
#o(Q) -#o(0) in (31) by its value at Q = e, - 
leading to 



■ m 



kin 



2E C 
3AEj 



Ej_ 
E c 



O 



Ej_ 
E C 



(33) 



where = ti 2 /Ej is the "naive" tight-binding mass 

of a single Cooper pair. The polaronic reduction of the 
mass is evident from (33) in the regime AEj > Ec > Ej. 



G. Adiabaticity condition and the validity of the 
mean field theory 

The analysis above rests on two assumptions: a) the 
dynamics of Q is slow and allows us to neglect the 
Landau-Zener tunneling in the derivation of Eq. (22); b) 
the field Q can be regarded as classical. 

Since in terms of Q the solitons are large objects (with 
the size of the order of 2A) one can expect the second 
assumption to hold in a wide parameter range. In par- 
ticular, at Ej ~ Ec we can estimate the "effective Planck 
constant" 20 for the Lagrangian (26) as 



/3 




< 1 • 



(34) 



While at small enough Ej / Ec strong suppression of the 
nonlinear Bloch inductance Lb(Q) for Q ^ e may be- 
come important, we expect this effect to be of minor 
significance in the intermidiate range of Ej /Ec- 

The situation with the assumption a) is much more 
tricky. First of all, the adiabaticity imposes an upper 
boundary for the velocity of solitons to be considered. In 
the limit Ej <C Ec the probability of a Landau-Zener 
transition into the first excited level is given by 



-A 2 J- 



(35) 



8 



where A « Ej, e(Q) is the difference of the charging 
energies of the two charge states involved in the process. 
Thus, i = QOqE and \8q£\ « If we demand P < 
e _:c , the corresponding limitation on the soliton's velocity 
reads 



2^2 



1 



2 Bc 



(36) 



We see that for Ej/Ec — >• 0, the maximal velocity for 
which adiabaticity still holds goes to zero. 

It is obvious that even for a static soliton the adia- 
baticity condition can be broken by fluctuations around 
the saddle point. Thus, the precise determination of the 
applicability region for the adiabatic approximation re- 
quires understanding of the characteristic time scale for 
the two-point correlation function of Q with the account 
of the nonlinear Bloch inductance. However, the good 
agreement between the mean-field theory and the tight- 
binding approach found in the calculation of the soliton 
mass for intermediate Ej/Eq ~ 0.5 allows us to expect 
that adiabaticity indeed holds in this parameter range 
for small soliton velocities. 



V. CONCLUSIONS 

In this paper we studied the dynamical properties 
of the charge carriers (charge solitons) in infinite one- 
dimensional Josephson arrays without disorder. We ap- 
plied two complementary techniques and arrived at our 
main result: in the parameter regime Ej < Ec < AEj 
the polaronic effects strongly reduce the effective mass of 
the charge solitons which scales approximately as EJ 2 . 
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Appendix A: Low lying excitations and the spin 
formulation 



The aim of the present Appendix is to find an explicit 
description of the low energy charge configurations satis- 
fying the constraint 



^ \r — r'\n r n r > = . 



(Al) 



Let us consider one such configuration. We try to add a 
dipole at islands R and R + 1 to this configuration, i.e. 



we construct a new configuration given by 

fl r = U r + S rR - 5 ryR +i . 



(A2) 



For the new configuration to belong to the low lying sec- 
tor we need the following condition to hold: 

\r-r'\n r h r > =2^n r {\r-R\ - \r-R- l|)-2 = . 

rr' r 

(A3) 



We can rewrite this condition as 

n r — n r = 1 . 

r>R+l r<R 



(A4) 



In terms of the spin variable cj z introduced in Sec. IIIB, 
Eq. (A4) is equivalent to the condition 



°r = "I • 



(A5) 



After the creation of an additional dipole (i.e. the Cooper 
pair tunneling from island R-\-l to island R) the new state 
is described by 



h r> + J2 



■2& 



r,R • 



(A6) 



r'>r+l 



We thus conclude that the Cooper pair tunneling from 
island R + 1 to island R is allowed (i.e. drives the sys- 
tem into another state within the low energy subspace) 
if o z R = —1. Such a tunneling corresponds to the spin 
flip at the link connecting the island R and R + 1. It is 
easy to check that the inverse process (tunneling from R 
to R + 1) is allowed only when a z R = 1 and also leads to 
the flip of (j z R . Taking into account that a single Cooper 
pair at island R is described in terms of a z by a domain 
wall 



J r<R 



-1 



} r<R 



1 , 



(A7) 



we conclude that the low energy configurations are those 
with cr z = ±1 for all r Translated into the charge lan- 
guage this condition gives the conditions mentioned in 
the main text (Sec. IIIB). 



Appendix B: Projecting the Hamiltonian 

We consider the action of the Josephson term in the 
Hamiltonian (8) on the state (9). Obviously, we can drop 
all the terms with r > R + wq or r < R — wq — 1 from the 
sum over r since acting on the state (9) they inevitably 
create the configuration of the width greater than wq. 
The terms with R + 1 < r < R + wq do not change 
the position of the first spin up in the chain and, thus, 
the coordinate of the soliton R. Thus, their action is 
described by 

R+wq w 



<|Wl,...^ ) = I^^Sflcri,...^) . 

(Bl) 



r=R+l 
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We consider now the action of cr^. It is convenient to 
introduce the operator of the right cyclic shift T acting 
on the states \d\ , . . . cr Wo ) according to 



T\cti , . . .a WQ ) = \a Wo , cti , . . . <j 2 ) • 



(B2) 



Assume that in the state \5\ , . . . <r Wo ) exactly /c first spins 
are down ( — 1) (we require now < k < wq — 1; the case 
of i^o spins down will be considered separately). The di- 
rection of other spins is arbitrary. In this case the action 
of cr^ on our state is given by 



Or\R)\&1,...0wo) 

= \R + k + l) [T+] 



fc+i 



fc+i 



cr fc+1 |cri ,-.-Cr«;o> • 

(B3) 



On the other hand, if in the given state exactly m ^ k 
first spins are —1, then 



fe+i 



a k+1 \ai,...a Wo ) =0 . 



(B4) 



Thus we conclude that for any state (except the state 
with all spins down) 



a x R \R)\a u ...a Wo ) 

Wq — 1 

= \R + k + l) [T + ] 

k=0 



fe+1 



Finally, taking into account that 
o%\R)\l,...l) = \R+wo + l) 



n 



we find for the projection of a x R onto the space spaned 
by the configurations of the width less or equal than wq 



a%\R)\ai,...a Wo ) 



k 



|a-i , . . . a-^ > . (B7) 



The last part of the Josephson Hamiltonian 

R-w-l 



(B8) 



r=R-l 



after the projection onto the subspace of interest, pro- 
duces an expression conjugate to (B7). Thus, summing 



up all the contributions we find the projected Josephson 
Hamiltonian 



k=l 



Hj\R)\a 1: ...a Wo ) = -Ej\R)^ j dl\d 1 ,...d WQ ) 

h 

w 

J2\R + k)[T + ]' 



E, 



fc=i 













w 






R + w + 1) 


ik 


+ /i.e. 


01 , • • •<) 











CT fe+1 |CTi,...CT Wo ) . 

(B5) 



I,... I) (B6) 



Going to the momentum domain with respect to the 
cyclic coordinate R and performing some simplifications 
based on the elementary properties of T and <r^, we fi- 
nally arrive at Eq. (13). 



Appendix C: Adiabatic calculation 

The most direct way is to apply the time-dependent 
perturbation theory. We will instead use a technique pro- 
posed by Berry 2 where we apply several unitary trans- 
formations, so that the eigenstates of the transformed 
Hamiltonians approach asymptotically the actual evolv- 
ing state. 

We consider the adiabatic case and the terms of orders 
higher than Q, Q 2 are omitted. The general transforma- 
tion reads 



H k+1 (t) = U J k H k U k -ihUlU k , 



(CI) 



with U k (t) being the time-dependent unitary operator 
that diagonalizes H k at each time t. We can write 



ei k Ht)) = U k (t)\n) , 



where effl (t) ) are the eigenvectors satisfying 



H k {t) e«(t)\ = £&*>(*) e«(t) 



(C2) 



(C3) 



and time-independent vectors \n) can be chosen arbitrar- 
ily, e.g., \n) 



oo)y. We now perform the first 
step of this process explicitly. As Hamiltonian H G (t) we 
take (21) with ei 0) (t)\ = \e n (Q(t))). From (C2) we find 



(C4) 



Applying (CI) gives the transformed Hamiltonian 



ffi(t) = 5>> (n\E^-ihJ2\n)(e^ 



m,n 



i(0) 



(m\ . 



(C5) 
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We find 



(t) j by applying the usual time-independent 

perturbation theory with the perturbation being the sec- 
ond term on the RHS of (C5). This gives the new trans- 



formation matrix U\(t) because 



e { n\t) 



Ui(t)\n), 



£n\t)^ we can go back to the desired en\t) 
basis via 



o(°)/ 



K) 



(C6) 



which can be used for a second transformation to ob- 
tain H2 (t) with eigenvectors effl (t) ) . After calculating Putting everything together gives 



!*(*)> = K) 



It 



2 



fmn(t) 

(*) 



km) 



fmkfkn frr 



Qm.n, ^ 



kn 



tt 2 



frnn d^l mr 



(C7) 



with f nm Q (d Q e n \ e m ) and Vt nm En h Ern . 
We are now able to calculate the voltage V r with 
where we set n = in (C7) as we consider our system be- 



ing initially in the ground state. We use (e n | 3qH \ e m ) = 
) for n ^ m and (e n \ d Q H \ e n ) = d Q E n 
and arrive at Eqs. (22,23). 
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